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Two antiferromagnetically coupled tunneling systems is a minimal model exhibiting the effect 
of quantum-mechanical phase in the Landau-Zener effect, ft is shown that the averaged staying 
probability oscillates vs resonance shift between the two particles, as well as vs sweeping rate. Such 
a resonance shift can be produced in Mn4 dimers by the gradient of the magnetic field. 
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Tunneling at an avoided level crossing, or the Landau- 
Zener (LZ) effeclpii^ is a quantum phenomenon that was 
much studied in physics of atomic and molecular colli- 
sions. Recently an experimental technique using the LZ 
effect was applied to single-molecule magnets to extract 
their tunneling level splitting A. 3 - 4 

In spite of the quantum nature of the LZ effect, its 
basic form can be described classically by a Landau- 
Lifshitz equation for a magnetic moment in a time- 
dependent fieldA^ However different kinds of interac- 
tions between tunneling magnetic molecules in a crystal 
make the LZ effect much more complicated. If the inter- 
actions are treated within the mean-field approximation 
(MFA), then the LZ tunneling can still be described by 
a nonlinear Schrodinger equation or, equivalently, by a 
classical nonlinear Landau-Lifshitz equation. In general, 
one is left with a tremendous problem of solving a full 
Schrodinger equation for an A-particle system. 

The question of how good is the MFA for the LZ ef- 
fect with interaction was studied in detail for the ide- 
alized "spin bag" model of N tunneling particles inter- 
acting each with each with the same coupling strength 
j | 6 4 7 This model can be mapped onto the problem of the 
giant spin S = N/2, so that the MFA limit N — > oo cor- 
responds to the classical limit S — > oo. In Ref. 7 it was 
shown that if the classical trajectory is smooth, the MFA 
yields qualitatively correct results, and quantum correc- 
tions can be calculated with the help of the cumulant 
expansion. In the case of complicated classical motion 
the MFA becomes unreliable. 

For the spin-bag model, both the MFA and the full 
quantum-mechanical solutions show that the ferromag- 
netic coupling suppresses LZ transitions (i.e., increases 
the LZ staying probability P), whereas the antiferromag- 
netic (AF) coupling increases transitions. This is in ac- 
cord with physical expectations based on the time depen- 
dence of the total field on a magnetic molecule, the sum 
of the external sweep field and the molecular field from 
the neighboring molecules. If one of the molecules tun- 
nels, then for the ferromagnetic coupling the jump of the 
total field is positive, so that the neigboring molecules are 
being brought past the resonance and lose their chance 
to tunnel. For the AF-coupling, the jump of the total 
field is negative, so that the neigboring molecules recieve 
an additional chance to tunnel and thus P decreases. 

A more realistic model for the LZ effect with in- 
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FIG. 1: Bare energy levels for two antiferromagnetically cou- 
pled tunneling systems for V < 4J. 



teraction should incorporate both distance-dependent 
couplings and individual resonance fields for magnetic 
molecules. Recently obtained solution for this model in 
the fast-sweep limits showed an oscillating dependence 
of the averaged one-particle staying probability P avr of 
the system on the resonance shifts between the molecules 
on the ith and jth lattice sites. This is an essentially 
quantum-mechanical effect arising due to the possibil- 
ity to reach the same final state by different sequences 
of individual tunneling events. The different quantum- 
mechanical phases accumulated on different ways lead to 
the interference in the final state. This poses a major 
challenge in the theoretical description of macroscopic 
systems since the phase factors are very sensitive to the 
microscopic details and dephasing processes should play 
a big role. Certainly the effect of interfering tunneling 
paths cannot be described by the MFA. 

The minimal model that exhibits the phase effect in the 
LZ tunneling is the model of two antiferromagnetically 
coupled tunneling systems (see Fig. that describes 
a particular transition in the recently discovered Mn-4 
dimeri2ii& We will use the Hamiltonian 
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H = ^,{[W{t) - VS] a iz + Av lx } + Ja lz <j 2z , (1) 
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where <Ti are Pauli matrices, W(t) — vt is the global 
time-linear energy sweep, A is the level splitting, J is 
the coupling, and we set 



-V/2, V 2 = V/2, 



so that V = V2 — Vt is the resonance shift. For V > 
particle 1 crosses the resonance first. The bare energy 
eigenvalues for this Hamiltonian are 
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~W{t) + 

V/2 - J, 



(2) 



where +, — means spin 1 up and spin 2 down. There are 
four first-order level crossings: 



W(t) = -V/2 - 2 J = Wt 
W{t) = V/2 -2J = W 2 
W{t) = -V/2 + 2J = W 3 
W{t) = V/2 + 2 J = W 4 . (3) 



The central crossing, £L._ = E +i+ =>■ W = is the 
second-order crossing with the splitting ~ A 2 / J that can 
be neglected in the case of well-separated resonances^ In 
this case one has four independent LZ transitions, each 
described by a scattering matrix (see, e.g., Ref. Illrt 
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where 



P = 



7rA 2 

2hv 



(4) 



(5) 



is the Landau-Zener staying probability and 

= 7r/4 + Argr(l-i<5) + 5(ln<5-l) (6) 

with 5 = ej (2ir). Evolution of the wave function between 
level crossings reduces to the accumulation of the phase 
factors exp [i$ a p(t)] , where the phases 



*a/j(t) 



dt'E al3 (t'), a,P = ± (7) 



can be easily calculated for the linear sweep from Eq. (0) 
and are quadratic in W. The wave function of the system 
can be written as ^(t) — YlapCafiit) \ a P) ■ Before the 
first crossing one has both spins down, C- y - = 1 (drop- 
ping an irrelevant phase factor) and C a p = otherwise. 
We use thus defined wave function C^g as the initial con- 
dition and we introduce 

^ a p{W 2 ,W 1 ) = ^ a p{W 2 )-^ a [3{W 1 ). (8) 

Then after the fourth crossing one has 

<W = (8a», +M f3 „ , + 5 „,_ 

xe iA^, l0 ,{W M (Sp' i+ M a >i a i + 5p, _5 a » a >) 



that is the final state. In this formula summation over 
repeated indices is implied. The staying probability for 
the first and second particles are given by 



Pi 
Pi 



(10) 



The average one-particle staying probability and the re- 
duced magnetization read 



P. 



avr - \{Pt+P 2 )=l-\\C 



+ ,- 
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M = l-2P avr . (11) 
After some algebra one obtains from Eq. © the results 
P 1 =P 2 {2-P) (12) 
P 2 = P(2-P)(l-P + P 2 )-2P(l-P) 2 cos 
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and 



P 



Pavr - P ( 1 - ~ ) (1 + P 2 ) - P(l - P) 2 COS 
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where the argument of the cos can be rewritten as 
UV 8JV 
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(13) 



(14) 



This is exactly the phase argument in Eq. (8) of Ref. 
8. The oscillating quantum-phase term in our solution 
arises because the state \H) can be reached in two dif- 
ferent ways: (i) Spin 1 flips first and spin 2 flips second 
(crossings at W = Wt and W4); (ii) Spin 2 flips first 
and spin 1 flips second (crossings at W = W 2 and W3). 
The amplitudes of these two processes add up, and the 
accumulated phase difference leads to oscillations. 

The interference effect in a system of two tunneling 
particles takes place for the AF coupling only. For the 
ferromagnetic coupling, the two horizontal lines corre- 
sponding to |t|) and ||t) g° above the (inactive) cental 
crossing, cf. Fig. Q As a result, there are only transi- 
tions at crossings at Wt and W 2 [see Eqs. (0] but no 
transitions at W3 and W4. Instead of Eq. © one has 

(15) 

and the results for the probabilities are 



Pi = P, P 2 = l-P + P 2 

Pavr = (l+P 2 )/2. 



(16) 



Note that P avr coincides with Eq. (20) of Ref. |6| for 
N = 2 and is independent on the resonance shift. In this 
model P avr > 1/2 because tunneling of both particles is 
impossible, C+,+ = 0. The case of a strong resonance 
shift corresponds to V > 4|J|. In this case the coupling 
plays no role, and one obtains P\= P 2 = P. 

Let us now consider the fast-sweep limit, e <ti 1. In this 
case one can write the expansion of the averaged staying 
probability in the form 



P, 
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72 + £ 2 /o, 



(17) 
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where Iq describes the deviation from the standard LZ 
effect, Eq. due to the interaction^ For not too large 
resonance shifts, V < 4| J|, one obtains from Eqs. 
and lUnj) 



h = 



-1/2- 
1/2, 



' sjv , 



J > 
J < 



(18) 



that is equivalent to Eq. (18) of Ref. |3- For large cou- 
plings and resonance shifts, the cos-term in Eq. H18|) 
oscillates fast and averages out. In this case one can 
conclude that the effects of antiferro- and ferromagnetic 
couplings are just the opposite: The former enhances 
transitions while the latter suppresses transitions by the 
same amount. 

Our present model allows, however, to analyze the in- 
fluence of interactions in the whole range of sweep rates, 
and it shows that the effect of the AF coupling is smaller 
than that of the ferromagnetic coupling. Dropping the 
cos-term in Eq. (|13|) , one obtains 
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P 



-P 2 (l 



P) 2 /2, 



(i-iy/2, 



J > 
J < 0. 



(19) 



For J > this difference is small everywhere, and its 
absolute value attains a maximum at P = 1/2, where 
Pavr — P = —1/32. On the other hand, for the ferromag- 
netic coupling P avr — P tends to 1/2 in the slow-sweep 
limit, £ ^> 1. 

Numerical solution of the LZ problem for two AF- 
coupled tunneling systems is shown in Fig. [2J One can 
see from Fig. that for J 3> A the staying probabilities 
P2 and P avr begin to oscillate starting from the values of 
the resonance shift V that satisfy V <C A. The numerical 
data for P avr is well described by Eq. (|T3|l starting from 
V ~ 0.1. On the other hand, for Pi and P2 the condi- 
tion of well-separated resonances is more restrictive and 
it requires somewhat greater values of V to validate Eqs. 
(|12fl . Fig. |3» shows oscillations of P2 as a function of 
the sweep parameter e, as well as a faster decay of Pi in 
comparison to the standard LZ effect. 

Let us now discuss the application of the present re- 
sults to Mn4 dimers. The coupling between the two 
Mn4 monomers with S = 9/2 was shown to have the 
form of the isotropic exchange^ H ex = J cx Si • S2 
with J cx ~ 0.1 K for the mostly studied compoundAifi 
Density-functional theory calculations 13 yield somewhat 
larger values of J ex . The uniaxial anisotropy D ~ 0.75 
K creates a barrier for spin tunneling. The level split- 
ting between the ground states | ±9/2) in the Mn4 
monomer is A ~ 2xl0~ 7 K (Ref. Q) and it should 
be of the same order of magnitude in the Mn4 dimer. 
The ground-state tunneling in a Mn4 dimer can be de- 
scribed by a pseudospin Hamiltonian of Eq. JIJ with 
J = S 2 J CX ~ 2 K. The period of oscillations on the res- 
onance shift V follows from Eqs. (|13JI and (|14|l and is 



given by Vp er iod = ( 7r2 /4) A 2 /(Je). For e — 1 one ob- 
tains Vpcriod — 5 x 10~ 14 K that corresponds to the dif- 
ference of longitudinal magnetic fields AH ~ 10~ 14 T be- 
tween the two monomers. With the distance between the 
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FIG. 2: Staying probabilities for a system of two antifer- 
romagnetically coupled particles o - vs. resonanse shift V. 
Asymptotic values are shown by horizontal dased lines, b - 
vs sweep parameter e. 



monomers I ~ 10 A, this amounts to the very small field 
gradient, dH/dx ~ 10~ 3 Gauss/cm! This means that 
small inhomogeneities of the magnetic field, due to, e.g., 
dipole-dipole interaction, should average out the quan- 
tum oscillations in the LZ effect. For possible applica- 
tions in quantum computing, it is desirable to have a 
larger tunnel splitting A, to achieve a faster performance 
rate and reduce the influence of decoherence. This can 
be achieved by applying a transverse magnetic field. In 
this case the period of the quantum oscillations consid- 
ered here will be much larger, and their observation will 
require much robuster field gradients that will exceed un- 
controlled inhomogeneities of the magnetic field. 

The author thanks R. Schilling for many stimulating 
discussions. 
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